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Abstract. We review the method of the calculation of multiloop integrals suggested in Ref. [T]. 



Qh' 1. Introduction 

pi |! The calculation of the multiloop integrals is important for the applications of the perturbative 

^ ' quantum field theory. A number of the calculational methods has been developed for this 

purpose. They roughly fall in to two categories: direct and indirect methods. The former 
include Feynman parameterization technique, Mellin-Barnes representation method, Gegenbauer 
%^J ■ polynomials expansion in coordinate space. The latter methods make use of the equations which 

QQ , are satisfied by the integrals. The construction of these equations relies on the possibility to 

\0 ] make the IBP reduction of the integrals under consideration to some finite set of the integrals, 

called "master integrals". Using this technique one can derive the differential equation for 
the integrals with several scales, or difference equation for the integrals with one scale. These 
equations determine the integrals up to the solution of their homogeneous parts. 



m 

^vq , In Ref. m a method of multiloop integrals evaluation based on D recurrence relations [2] 

and P-analyticity was suggested (DRA method). Since then it was applied to the calculation of 
various complicated integrals O IU O El [71 [8] . In this contribution we review this method. 

*^ ; 2. DRA method 

.^.' 2.1. Dimensional recurrence relation 

Assume that we are interested in the calculation of the L-loop integral depending on E linearly 
independent external momenta pi, . . . ,pE. There are A^ = L{L + l)/2 + LE scalar products 
depending on the loop momenta li: 



^l] — oji 



h-Qj-, i = l,...,L; j = l,...,K, (1) 



where gi,...,L = h,...,L, qL+i,...,K = Pi,...,E, and K = L + E. 
The loop integral has the form 

J^""^ ^^^- ■ -^^) = j :;^ij4^t^^^^ (2) 

where the scalar functions Da are linear polynomials with respect to Sij. The functions D^ are 
assumed to be linearly independent and to form a complete basis in the sense that any non-zero 



linear combination of them depends on the loop momenta, and any Sij can be expressed in terms 
of £>„. 

The integral J^ ' (ni, . . . , n^v) = J^ ' (n) can be considered as a value of function of integer 
A^-dimensional vector n = (ni, . . . jUn). The set of points in Z (and the corresponding set 
of integrals), having the same set of positive coordinates, form a sector. Thus the whole space 
is split into 2 different sectors which can be conveniently labeled by their simplest elements. 
E.g., (1, 1,0,0, .. . ,0) denotes the sector in which all points have exactly two first coordinates 
positive. The points (and the corresponding integrals) in Z^ and the sectors form a partially 
ordered set with respect to the following relation. We say that the point n is simpler than the 
point m iff all coordinates, positive in n, are also positive in m, but not vice versa. 

It is convenient [T] to introduce the operators Ai,Bi which act on such functions as 

{Aif ){..., Hi,...) = nif{...,7ii + l,...), 

(i?,/)(..., n„... ) = /(..., n^-l,...). (3) 

These operators respect the notion of sectors in a sense that the values of functions Aif and 
Bif in the point of some sector are expressed via the values of function / in the same or simpler 
sectors. Using these operators, we can express the integral in P it 2 dimensions via the integrals 
in P dimensions [9]: 



(V-E-L + l)^ 

j(^-2)(n) = (^/2)^Q(Ai,...,A^)j(^)(n) , (5) 

where /i = ±1 for the Euclidean/pseudoEuclidean case, and polynomials P and Q are determined 



as 



P{Di,...,Dn) =det {qi-qj} 



i,j=l...K ' 



These relations are the most useful when the integral in the left-hand side is master. Then, 
making the IBP reduction of the right-hand side of Eqs. dH,©, one obtains difference equations 
for the master integrals. After this, we obtain a linear combination of the integrals of the same, 
or simpler, sectors as the integral in the left-hand side of Eqs. (I1D,([5]). We assume that the 
integrals of the the simpler sectors are already known at this stage, either by the same, or by 
some other method. 

The IBP reduction makes Eqs. ([I]),© equivalent, so in what follows we consider only the 
lowering dimensional recurrence relation, obtained from Eq. ^. We also change notations, 
writing 

J(^) = J(.), (7) 

where z^ = T>/2. If a sector of the integral in the left-hand side contains only one master integral, 
the general form of the equation is the following: 

J{„ + l) = C{i^)J{i^) + Rii^), (8) 

where C{i') is some rational function of v and R{iy) is a non-homogeneous part constructed of 
the simpler master integrals in T> = 2]^ dimensions. If a sector contains several master integrals, 



it is convenient to consider them as a column- vector integral J(z^) and to write the equations 
for them in a matrix form: 

J(zy + l) = C(i/)JH + R(zy), (9) 

where C(z^) is some rational matrix (i.e., the matrix with coefhcients being rational functions of 
v). 

2.2. General solution 

In order to write the general solution of this equation, we introduce the notion of indefinite sum 
and indefinite product. Given the function F{v), the indefinite sum TjF{v) is any function, 
satisfying the relation 

^F{u + l)-T.F{v) = F{u) (10) 

The indefinite sum S F[v) is defined up to the addition of arbitrary periodic function of u. If 
the function F{v) decreases faster than 1/v when v — t- +oo and/or u — t- —do, one can write 



S F{v) = S F{v) + uj{z) = -YF{u + n)+ uj{z) (11) 



def 

n=0 

and/or 

-1 

^F{u)= ^ F{u)+uj{z)= Y F{v + n)+oj{z), (12) 



where oj{z) = ti;(exp(27riz^)) is arbitrary periodic function. 

The indefinite product I[F[v) is any function, satisfying the relation 

IlF{v + l)/\lF{v) = F{v) (13) 

YiFiy) is defined up to the multiplication by arbitrary periodic function of u. Note that for 
any rational function F{i') one can express YiF{v) as a product of F-functions multiplied by 
arbitrary periodic function, see Ref. [1]. We will determine the indefinite product nF(i/) also 
for the matrix functions F(i/) as a matrix function, satisfying the relation 

IW{v + l){YW{v))-^ =¥{v) (14) 

Obviously, this defines nF(z/) up to the multiplication from the right by arbitrary periodic 
matrix. 

In these notations one can write down the general solution of Eq. ([9]) in the form 

3{v) = E{v)-^^E{v)'K{v), (15) 

where S(i^) = [nC(z/)]~^ is called a summing factor. If the summand in Eq. (|16p decreases faster 
than 1/u when v — >• +00 and/or v — >• —00, one can use Eqs. pT]) . (fT2|) to write 

J(z^) =S(z^)"^cj(z) + S(i/)-i S E{v)'K{u). (16) 

±00 

Note that the second term in the right-hand side is invariant with respect to the choice of the 
summing factor. Thus, in order to find the integral, we need to fix the function cj(z) for any 
choice of the summing factor. Here we note the essential difference between the general solutions 
of differential and difference equations. The homogeneous part of the general solution of the 
differential equation is parameterized by finite number of constants, while the homogeneous part 
of the general solution of the difference equation is parameterized by the periodic function(s). 



■J'^0 + 1 




Figure 1. Mapping v ^ z 

While fixing the finite number of constants is a relatively simple task, fixing the periodic function 
may be quite difficult. The key point of the DRA approach is the usage of the analytical 
properties of the integral in order to fix this function. Let us rewrite Eq. (|16|) as 



u{z) = S(z^)J(i/) 



S S(i/)R(z>). 

±00 



(17) 



Suppose that we can show that the right-hand side of this equation is a meromorphic function 
of z in the Riemann sphere. Then ijj{z) is determined, up to a constant, by the principal parts 
of its Laurent series around singular points. The analysis of the analytical properties of the 
right-hand side of Eq. ()17p can be performed in any vertical stripe of unit length in the complex 
plane of v, as demonstrated in Fig. 1. There is no cut on the ray arg(z) = 2t:uq because the 
right-hand side is necessarily a periodic function of v. The proper choice of this stripe {basic 
stripe) can simplify the evaluation of the integral. 



2.3. Analytical properties of integrals 

Usually, one is interested in the analytical properties of the integral as a function of some external 
invariant. Here we are interested instead in the analytical properties of the integral as a function 
of z/ = D/2. It appears that the latter properties are much simpler than the former. 

Let us consider the parametric representation of some L-loop integral in Euclidean space with 
/ internal lines (see, e.g., Ref. [TO]): 



J {v) = r (/ — Lv) I dxi . . . dxi6 (1 — Yli^i) 



[q (x)] 



vL~I 



\p{x. 



^u{L+l)~I■ 



(18) 



The polynomials q {x) and p (x) are determined by the graph. For our consideration it is 
important only that both these functions are nonnegative in the whole integration region, 

q{x)^Q, p{x)^0. 

Suppose that the parametric representation converges for all v in some interval (z^i,^'2)- Then 
it is easy to see that it converges on the whole stripe S = {u, Rez/ E (z^i, 1^2)}- Indeed, we can 
estimate 

Therefore, Eq. ()18p determines J(z^) as a holomorphic function on the whole stripe S. In fact, 
in Ref. |11] a more general fact was proved: the representation (jlSp determines the meromorphic 







n., (c.l = 0) T4,i (c.l. = 1) r4,2(c.L = 1) n.^ (c.l. = 2) Tfij (c.l. = 3) 
Figure 2. Three-loop all-massive tadpole master ntegrals. 

function on the whole complex z^-plane. The determination of poles position can be automatized 
thanks to FIESTA [T2] . 

As Im I' tends to ±cx) while Re u is fixed, Eq. (|19p shows that 

\J{u)\ < const X e"^^^ |Imi/|^"^/2-LRc^ ^ ^^^^^ ^ |z|^(^+°)/^ . (20) 

Thus, the parametric representation can be used for the determination of the analytical 
properties of the first term in the right-hand side of Eq. p7|) . Since the second term in the 
right-hand side of Eq. p7|) is constructed of the integrals already known at this stage, we, 
indeed, can determine the position of the poles of function co{z). 

2.4- Numerical issues 

In order to demonstrate how the representations obtained within the DRA method can be 
treated numerically, let us consider the three-loop master integrals depicted in Fig. [2l Direct 
application of the DRA method gives the following results: 

n,i = r(i - i^f (21) 

^''' = i:A^~-''o^^'^'^'''}' ^''' ^ r(2 - 2' )r(i - u) ('') 

J. _ 1 /3vr2 ^ (lli^ - 8)S4,2T3,i ) 4i-3.r(2 - z.)^ 

'''"54:^11^"-- W^ /' '''"r(3-3z.)r(3-2i.) ^^'^^ 

1 ri67r2 4(7j, _ 4)T3,^ - 36(3z. - 2)r4,i - (2 + 13i/)T4,2 

-tS,! — t; — S —^ r Zj 05,1 



S5,i i 27 -^ "' 27(31/ -2) 

9-'^r(l- ;y) 
r(2-2i/ 



55,i(^) = ^.. ' ...2 (24) 



2 , ^ 'S'6,1 /m om n1 rt 2 "^ 



The above results for the integrals T^^i and Tg^i contain repeated sums of the form: 

00 00 00 

^=EE--- E /i(™i)/2M.../fc(nfc) (26) 

ni=0n2=ni nfc=nfc_i 

From the viewpoint of numerical calculation, it is important that the dependence on the 
summation variables is factorized in the summand. We note that it is not the case for the 
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Figure 3. Time of calculation (sec) as a function of (a) requested e-order (T6^i(4 — 2e) with 
200-digit precision); (b) requested precision (r6^i(4 — 2e) up to 0(e'^)); (c) complexity level 
(expansion up to O(e^) with 500-digits precision) 



Mellin-Barnes multiple sums. Let us show that the calculation of such factorized sums can be 
performed without nested loops. Indeed, let us define k sequences si(n), . . . , Sfc(n) as follows: 

Si{-1) = 0, si(n ^ 0) = si(n- 1) +/i(n- 1), Sj>i(n ^ 0) = Si{n - 1) + fi{n)si-i{n) 

(27) 

It is easy to check that the last sequence, Sk{n), converges to s. Obviously, the iterative 
calculation of these quantities can be organized in one loop. Therefore, the complexity of 
the numerical calculation grows moderately with the complexity level of the integrals. The 
dependence of the calculation time on the order of the expansion, required precision, and 
complexity level of the integrals is shown in Fig. [3j We see that for all three cases the method 
shows quite moderate growth of the evaluation time. 

3. Conclusion 

We have briefly reviewed the method of calculation of multiloop integrals based on the P- 
recurrence and P-analyticity. The method appears to be powerful enough to deal with the most 
complicated cases. 
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